Search results for "numerical [Methods]"
showing 10 items of 500 documents
A model study of Hartree-Fock and Linear Response in coordinate space
1979
A fast procedure for spherical Hartree-Fock is obtained by coordinate space representation and a modification of gradient iteration. Along similar lines, the corresponding Linear Response equations are derived and solved, in order to achieve a fully consistent treatment. The Linear Response equations are applied to a change in particle numbers, i.e. to the description of isotopic differences. In a model study we look for their physical and numerical properties, i.e. linearity of the response, numerical stability and consistency requirements for the Hartree-Fock basis.
Linear response strength functions with iterative Arnoldi diagonalization
2009
We report on an implementation of a new method to calculate RPA strength functions with iterative non-hermitian Arnoldi diagonalization method, which does not explicitly calculate and store the RPA matrix. We discuss the treatment of spurious modes, numerical stability, and how the method scales as the used model space is enlarged. We perform the particle-hole RPA benchmark calculations for double magic nucleus 132Sn and compare the resulting electromagnetic strength functions against those obtained within the standard RPA.
Combining spectral and shock-capturing methods: A new numerical approach for 3D relativistic core collapse simulations
2005
We present a new three-dimensional general relativistic hydrodynamics code which is intended for simulations of stellar core collapse to a neutron star, as well as pulsations and instabilities of rotating relativistic stars. Contrary to the common approach followed in most existing three-dimensional numerical relativity codes which are based in Cartesian coordinates, in this code both the metric and the hydrodynamics equations are formulated and solved numerically using spherical polar coordinates. A distinctive feature of this new code is the combination of two types of accurate numerical schemes specifically designed to solve each system of equations. More precisely, the code uses spectra…
Cross-phase modulational instability induced by Raman scattering in highly birefringent fiber
2013
We report experimental and theoretical studies of Raman-induced cross-phase modulational instabilities (XPMI) in a high-birefringence, normally dispersive optical fiber. Experimental results reveal that the Raman-Stokes wave, generated by a quasi-CW pump beam, interacts with the latter to create a novel type of XPMI sidebands. These sidebands are characterized by a narrow gain bandwidth. The sideband frequencies are well reproduced by a linear stability analysis as well as by full numerical solutions of the coupled generalized nonlinear Schrödinger equations.
Turbulence structure and budgets in curved pipes
2013
Abstract Turbulent flow in curved pipes was investigated by Direct Numerical Simulation. Three curvatures δ (pipe radius a /curvature radius c ) were examined: δ = 0 (straight pipe), simulated for validation and comparison purposes; δ = 0.1; and δ = 0.3. The friction velocity Reynolds number (based on the pipe radius a ) was 500 in all cases, yielding bulk Reynolds numbers of ∼17,000, ∼15,000 and ∼12,000 for δ = 0, 0.1 and 0.3, respectively. The computational domain was ten pipe radii in length and was resolved by up to 20 × 10 6 hexahedral finite volumes. The time step was chosen equal to a wall time unit; 1 Large Eddy TurnOver Time (LETOT) was thus resolved by 500 time steps and simul…
Noise delayed decay of unstable states: theory versus numerical simulations
2004
We study the noise delayed decay of unstable nonequilibrium states in nonlinear dynamical systems within the framework of the overdamped Brownian motion model. We give the exact expressions for the decay times of unstable states for polynomial potential profiles and obtain nonmonotonic behavior of the decay times as a function of the noise intensity for the unstable nonequilibrium states. The analytical results are compared with numerical simulations.
Scale-free relaxation of a wave packet in a quantum well with power-law tails
2013
We propose a setup for which a power-law decay is predicted to be observable for generic and realistic conditions. The system we study is very simple: A quantum wave packet initially prepared in a potential well with (i) tails asymptotically decaying like ~ x^{-2} and (ii) an eigenvalues spectrum that shows a continuous part attached to the ground or equilibrium state. We analytically derive the asymptotic decay law from the spectral properties for generic, confined initial states. Our findings are supported by realistic numerical simulations for state-of-the-art expansion experiments with cold atoms.
Dynamical formation of a hairy black hole in a cavity from the decay of unstable solitons
2016
Recent numerical relativity simulations within the Einstein--Maxwell--(charged-)Klein-Gordon (EMcKG) system have shown that the non-linear evolution of a superradiantly unstable Reissner-Nordstr\"om black hole (BH) enclosed in a cavity, leads to the formation of a BH with scalar hair. Perturbative evidence for the stability of such hairy BHs has been independently established, confirming they are the true endpoints of the superradiant instability. The same EMcKG system admits also charged scalar soliton-type solutions, which can be either stable or unstable. Using numerical relativity techniques, we provide evidence that the time evolution of some of these $\textit{unstable}$ solitons leads…
Simulations of Array Configurations for the Square Kilometre Array (SKA)
2010
The Square Kilometre Array (SKA) is a new generation radio telescope for the next decades, working at metre to centimetre wavelengths. The SKA will be operational at the same time than other new optical, X-ray and Gamma-ray telescopes. It is of extreme importance that the SKA becomes competitive and complementary to those instruments. An extensive study of technologies and possible configurations involved is needed to ensure the SKA will reach the design specifications. To compare imaging capabilities between different SKA configurations or between the SKA and other instruments, we have implemented figures of merit based on several characteristics of these instruments. In this work we are p…
Diffusion in active magnetic colloids
2013
Abstract Properties of active colloids of circle swimmers are reviewed. As a particular example of active magnetic colloids the magnetotactic bacteria under the action of a rotating magnetic field is considered. The relation for a diffusion coefficient due to the random switching of the direction of rotation of their rotary motors is derived on the basis of the master equation. The obtained relation is confirmed by the direct numerical simulation of random trajectory of a magnetotactic bacterium under the action of the Poisson type internal noise due to the random switching of rotary motors. The results obtained are in qualitative and quantitative agreement with the available experimental r…